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solution  of  the  equilibrium  equations  for  a  viewport  in 
a  deep  submersible  are  defined..   A  likely  solution  method 
is  chosen  and  investigated  in  order  to  determine  its 
applicability  and  to  learn  more  of  the  behavior  and 
inter-relationships  of  the  governing  equations.   The 
difficulties  encountered  in  the  chosen  solution  method 
are  explained  and  another  approach  is  suggested. 
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INTRODUCTION 

Since  Piccard's  Trieste,   almost  all  deep 
subraersibles  have  utilized  small  plexiglass  viewports 
for  observation  and  control.   The  governing  structural 
mechanics  equations  for  the  viewports  have  not  been 
solved  and,  therefore,  the  stress  state  and  predicted 
deflections  have  not  been  known.   The  interaction  with 
the  pressure  hull  and  its  reinforcement  has  also  mainly 
been  known  from  finite  element  analyses. 

The  principal  investigations  of  the  behavior  of 
plexiglass  viewports  were  experimental  tests  reported  by 
Stachiw.  '  '  '    These  tests  included  short  term  loading 
of  conical,  flat,  and  spherical  windows  as  well  as  long 
term  creep  of  conical  windows.   The  data  from  these  tests 
provides  sufficient  information  for  safe,  large  safety 
factor  design.   It  does  not,  however,  provide  insight 
into  the  stresses  in  the  window  support  or  detailed 
information  of  expected  viewport  performance  and  useful 
life. 

Since  plexiglass  is  inexpensive,  the  viewports  have 

been  greatly  over-designed  and  the  present  state  of  the 

art  has  been  adequate.   However,  for  future,  more 

sophisticated  vehicles  with  larger  viewports,  a  better 
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understanding  of  the  stresses  and  deflections  in  the 
viewports  and  the  resultant  interaction  with  the  pressure 
hull  will  permit  a  safer  and  less  expensive  design  and 
greater  insight  into  critical  design  aspects. 


BACKGROUND 

Most  viewports  presently  utilized  in  undersea  vehicles 
are  truncated  cones  mounted  in  heavily  reinforced  openings 
in  the  pressure  hull.   They  are  held  in  position  by 
mounting  rings  but  the  primary  seal  is  simply  from  surface 
to  surface  contact.   The  contact  surfaces  are  usually 
polished  to  ensure  a  good  seal  and  a  90  degree  cone  angle 
is  utilized.   This  angle  was  determined  by  Piccard  to  be 
the  best  design  and  has  been  utilized  successfully  by 
designers  since  then.   The  tests  by  Stachiw  have  supported 
the  opinion  that  the  90°  cone  angle  is  best  for  acrylic 
viewports . 

It  can  be  seen  that  the  viewport  is  similar  to  a 

thick,  flat,  circular  plate  with  slanted  ends.   The  stress 

distribution  is  axi-symmetric  and  therefore  circular 

coordinates  should  be  utilized  in  forming  the  governing 

equations.   For  the  coordinate  system  in  Figure  1,  the 

deformation  is  symmetrical  around  the  z  axis  and  therefore 

all  derivatives  with  respect  to  6  are  zero  and  shear 

stresses  t  „  and  xQ   are  zero. 
ru      b  z 

Timoshenko  and  Goodier   show  that  the  equations  of 
equilibrium  therefore  take  the  form: 


Figure    1 


4////tfni 


6a    <5t     'a  -aQ 

r  +   «  +  _JL_i  =  Q  (1) 

6r     6z      r 


6x     6a    t  „ 

rz  +  5.  +  _^E  =  o 

6r     6z    r 


At  present  most  submersibles  have  moderate  depth 
capabilities  and  therefore  do  not  need  excessively  thick 
viewports.   The  desire  for  good  visibility  therefore  results 

in  a  design  with  a  thickness  to  diameter  ratio  of  usually 

2 
less  than  .5  for  the  viewport.   The  fact  that  Stachiw 

only  tested  viewports  with  thickness  to  inner  diameter 

ratios  of  .631  or  less  illustrates  that  this  is  presently 

the  major  area  of  concern.   Consequently,  the  compatibility 

equations  can  therefore  be  simplifed  by  assuming  plane 

stress.   Although  this  is  not  completely  accurate,  it 

should  allow  an  accurate  prediction  of  the  behavior  of 

almost  all  desirable  viewport  designs. 

The  compatibility  relations  therefore  take  the  form: 

2       2  1    62(ar+a0  +az) 

V^a £  (a  -aQ)  +  -4-  - — 2 L_  =  0 

r    r2  v  r   8'    1+v        „  2 

or 

2       2  11  6(ar+a0+az) 

V  aQ   +  -4-  (a  -a0)  +  ^-  -  ^— —   =  0  (2) 

0      2  v  r   0     1+v  r      6r 
r 

2      1   62(ar+a0+Oz) 


.2       1  T    +   1   ■g2(qr+g6+CTz)    n 
V  Trz  "  T"  Trz  +  1+V     6r~Si         ° 


The  strain  components,  therefore,  become: 


Su         u  Sw 

£r  "  6r  ,   £0  ":  r   '      z  '  <5z 


6u    <5_w 

'rz  ~  6z    6r 


Equations  (1)  and  (2)  can  be  satisfied  by  expressing 
the  stresses  in  terms  of  a  stress  function  cf> .   The  stresses, 
therefore,  become: 


6  ,  n2A    <52<J\ 
or 


6   /,,r72A    1  $<)>. 


Oz 


<S  2  &2<S> 

Trz   =    6r    [d-v)V  U  ~  — f] 

6z 


subject  to  the  condition  that: 


V2V2*  =  0  (3) 


These  equations  apply  to  a  viewport  in  an  undersea 
structure  or  pressure  facility  provided  the  appropriate 
boundary  conditions  are  met.   The  first  and  most  obvious 
boundary  condition  is  that  the  normal  stress  on  the 
pressurized  face  is  equal  to  and  opposing  the  external 
pressure.   Obviously,  the  condition  is  reversed  on  the 
unpressurized  face  where  the  normal  stress  is  zero.   The 
shear  stress  on  the  flat  faces  is  also  zero  since  they 
are  free. 

The  boundary  conditions  on  the  slanted  ends  are  less 
obvious.   In  reality  there  will  be  a  contraction  and 
rotation  of  the  mounting  reinforcement  with  pressure. 
However,  since  the  reinforcement  is  very  thick  steel 
with  a  high  modulus  of  elasticity,  it  can  be  safely 
considered  as  being  rigid  compared  to  the  plexiglass. 
Since  the  viewport  is  not  physically  attached  to  its 
mounting,  it  is  free  to  slide  along  the  mounting  opposed 
only  by  friction  and  the  shape  of  the  mounting. 

Two  boundaries  on  the  behavior  of  the  viewport  at 
its  mounting  can  be  specified.   The  viewport  can  be 
considered  to  be  rigidly  fixed  to  its  mounting  or  it  can 
be  allowed  to  slide  freely.   In  the  first  case,  it  seems 
that  the  stresses  produced  will  be  higher  than  the  actual 
situation  and  the  displacements  will  be  lower.   In  the 
second  case,  it  seems  that  the  stresses  produced  will  be 
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slightly  lower  than  the  actual  situation  and  the  displace- 
ments will  be  higher. 

The  first  case  is  satisfied  if  all  displacements 
are  zero  at  the  boundary,  i.e. 

D0     H 

u=w=0atr=-y-  (J  +  z)  tan  9 

The  second  case  is  satisfied  if  the  displacement  normal 
to  the  mounting  is  zero  at  the  boundary,  i.e. 

2     2 


0     H 
u  cos  8  +  w  sin  0  =  0  at  r  =  —  -  (—  +  z)  tan  6 


In  seeking  solutions  of  equations  (1) ,  (2) ,  and  (3) , 
Timoshenko  and  Goodier  considered  a  polynomial  form  of 
the  stress  function  $   which  also  satisfied  the  equation: 


6R2    *^    R2Ctn   ^  +  ?6*2-° 


This  resulted  in  particular  solutions  of  the  form: 


n  _  n(n-l)   R2  „(n-2)  ,n(n-l)  (n-2)  (n-3)„4  „(n-4) 
'n  "  a0Lz     2(2n-l)  R   Z      +  2-4(2n-l) (2n-3)R   Z 


and   additional  solutions  of  equations  (1)  and  (2)  found 

2 
by  multiplying  the  particular  solution  by  R  where 

2     2     2 
R   =  r   +  z  .   Consequently,  the  following  possible 


solutions  can  be  found: 


*0  =  ao 


*1  =  alZ 

<j>2   =  a2[2z2   -  r2] 

<f>3   =.  a3[2z3   -  3r2z] 

4>4   =  a4[8z4    -  24r2z2    +    3r4] 

4>      =  ac[8z5   -  40r2z3   +    15r4z] 

A                 r16    6  ..    4    2    ,     _A    2    4         5    6, 

q>      =  a^  [-=-z      -    40z   r      +    30z    r      -   Tr    ] 

6            6      o  -5 


and 


<J>2   =   b2[r2   +   z2] 


<J>3   =   b3[z3   +    zr2] 


4  2    2  4 

<|>4   =   b4[2z^    +   r"V    -   r*] 


5     2  3      4 
<j>5  =  b5[2zD  -  r  z   -  3r  z] 


4,  =  b^[8z6  -  16z4r2  -  21z2r4  +  3r6] 
6     o 


Timoshenko  and  Goodier  considered  a  flat  circular 
plate  loaded  with  a  uniform  pressure  on  one  face.  The 
boundary  conditions  were: 
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t    =  0  at  z  =  ±§ 
rz  2 


a      =  0  at  z  =  ? 
z  2 


a   =  -P  at  z  =  -J 
z  2 


/ 


These  are  identical  to  all  but  the  last  boundary  condition 
for  a  viewport. 

By  assuming  a  solution  form  of: 


A      ,16  6    An    4  2    ,«  2  4    5  6, 
<P  =  a,  (-^-z   -  40z  r   +  30z  r   -  -=r  ) 

DO  J 


+  br(8z6  -  16z4r2  -  21z2r4  +  3r6) 
b 


4       2  2      4 
+  a4(8z   -  24r  z'  +  3r  ) 


+  a3(2z3  -  3r2z) 


+  b3(z3  +  zr2) 


they  found  that 


0     3520HJ 


b   =  (— i^)  P 

352HJ 


11 


3 

a4  :   ("  384H}  P 


a    =  r1+5  V  )  P 

3    v72+60v;  r 


b3    ^~  24+20v)  P 


Therefore  the  stresses  become: 


or-PK2*v>£-  liftvi'A.s., 


3     -a       , 

_  pr    z   .££_£■] 

Z  orr      2  H     2 


x    =  3pr  (5-  -  z2) 
rz         4 


Using  the  equations  for  strains  given  previously  we 
find  that: 

1+v   6  2<{> 
u  =  ■  - 


E    6r6z 


w  =  ^  [2(l-v)^|  +  2(l-v)  |  ||  +  (1-2V)  i!|] 

<5r  6z 


Consequently, 


l+vr/2-v.  3  .  ,-3+3v.   3, ,-3.   ,.       5v   v  .  _ 
U  =  -j-I  (-3-)  z  r+  (—3-)  zr  +  (^)  zr+  (j^^)  r]  P 
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l+vr,-l-v,  4^,    3v.  2  2^,3-3v>  4  ,  ,  3  .  2..        3.  2 
W  =  "^-[(-^3)Z  +(^3)Z  r   (^3)r   W2   ("  8H)r 


+  (l2TlWzl  P 


The  final  boundary  condition  for  the  first  case  for 

D0     H 
the  viewport  requires  that  u=w=0atr=  -~-  -{■=■+    z)tan 

The  flat  plate  solution,  however,  gives 


1+V  _  r  ,-2+v,  4  .    fl  .  /3-3V,  4.3, 
u  =  -=—  P  [  ( =— )  z   tan  9  +  ( =■)  z   tan  ( 

E       H3  4H3 


,  ,2-v.  .3  ,  ,-9+9y.    3  ,   2, 

+  ( — ^)   A  z   +  ( ^— )  A  z   tan  f 

2H  8H 


,  ,9-9vv  A2  2  ,    n   ,  ,  3  ,   2  , 
+  ( 5-)  A  z   tan  0  +  (-rrr)  z   tan 

16H3  4H 


,-3+3 v*     .3  ,-3.     .       ,     ,      -5v^ 

~3~2~^~  Z  M  (12+W    Z    tan 


+  ^24+T5v)  A]  =  ° 


1+V  „  r  ,-1-v.   4    ,  3v,   4  .   2, 
W  =  -=-  P  [  ( — =— )  z   +  ( — =-)  7.      tan  f 

2HJ  2HJ 


+  <^)  z4  tan40  +  (=4)  Az3  tan 
16HJ  2H3 


13 


+  (Z3±3V}  Az3  tan30  +  (  3  ,  z2  +  (_3v  )    ^2 
8H3  8HJ 


+  (i|)z2  tan29  +  (i=^)  A2z2  tan2( 
8H  32HJ 


+   felo^    z  +   <M>   Az  tan  9  +   {^)   A'z  tan 


32H  256H3 


The  final  boundary  condition  for  the  second  case  for 

the  viewport  requires  that  u ' +  w  tan  0  =  0  at 

D0     H 
r  =  [ —  -  (-  +  z)  tan  6 ] .   The  flat  plate  solution, 

however,  gives 


u  +  w  tan  9  =  ^  P  [  (-^)  tan  6  z4  +  (^34)  tan39  z4 
E      2H2  4H3 


,3-3v,  .   5Q   4    ,2-v>  .3  ,   ,-9-3v.  ..   2a   3 
+  ( =■)  tan  6  z   +  ( — zr)    Az   +  ( — )  Atan  6  z 

16H  2H  8H 


+  (  3+3v)  Atan4G  z3  +  (^-34)  A2tan  0  z2  +  (~)    tan30  z2 
8H3  16H3  8H 
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+    (4)    tanez2+(^)A2tan30z2+(^±^)A32+(^±^)A3tan26z 
2il  32H  32H  32H 


+  ("  M)Az  +  (^)Atan2ez  +  (lrilo-)  tan9  z 


+    (^74)A4tan9+     (-3§h)     A2tan6+    <24W)A] 
256H 


where      A   =    (D      -   Htan  8  ) . 


Consequently,  it  can  be  seen  that  if  these  boundary 
conditions  could  also  be  satisfied,  solutions  would  result 
which  would  be  dependent  on  the  material  properties  v  and 
E  as  well  as  the  dimensions  Dn,  H,  and  0  and  the  pressure 
P.   Since  plexiglass  is  a  viscoelastic  material  the 
properties  v  and  E  will  be  slightly  non-linear  and  will 
be  functions  of  time.   In  practice,  all  viewports  are 
designed  to  remain  well  within  the  elastic  range  since 
vision  would  be  greatly  distorted  even  before  reaching 
the  non-linear  response.   Consequently,  a  solution 
assuming  that  the  material  is  linearly  elastic  will  be 
adequate  since  the  variation  in  properties  of  different 
batches  of  plexiglass  makes  more  exact  calculations  of 
little  use. 


SOLUTION  WITH  SIXTH  ORDER 
POLYNOMIAL  ASSUMPTION 

In  order  to  satisfy  either  of  these  additional 
boundary  conditions  it  is  apparent  that  supplemental  sixth 
order  terms  will  be  required.   Since  no  other  terms  are 
possible  using  the  method  of  Timoshenko  and  Goodier, 
other  forms  must  be  found.   The  desired  forms  must 

satisfy  equation  (3)  and  two  possibilities  therefore 

2  2  2 

become  apparent.   Either  VV<})  =  0orV<j>  =  0  will 

satisfy  equations  (1)  and  (2) . 

Probable  desired  terms  include  sixth,  fourth,  and 

third  order  terms.   General  terms  can  be  constructed  as 

follows : 

9  =  az   +  bz  r   +  cz  r   +  dr 

a      4,-22,    4 

9  =  ez   +  fz  r   +  gr 

a      3       2 
9  =  pz   +  qzr 

2  2 

In  order  to  satisfy  V  V  9  =  0  we  find: 

12b  +  8c 

a  ~       45 


,      3b  +  8c 
72 
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8      2  . 

e="3g  "  3f 
/ 

Additional  terms  can  be  constructed  as  follows 


<J>  =  hz6  +  iz4r2  +  jz2r4  +  kr6 


4  2    2  4 

<t>   =   lz      +  mz    r      +   nr 


2  2 

<J>   =   sz      +   tr 


2 
In  order  to  satisfy  V  <p  =  0  we  find: 


V,  16  V 

h  =  -  -=-  k 


i  =  24k 


j  =  -18k 


m  =  -8n 


i    8 
1  =  Tn 


The  boundary  conditions  for  the  combination  of  the 
flat  plate  solution  and  the  additional  terms  therefore 
become: 


-r=0atz=±§ 
rz  2 


a   =  0  at  z  =  J 
z  2 
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a     =   -P   at    z  =  -  t 
z  2 


u  =   w=Ooru  +   w   tan  9  =    0    at   r  "=   j^o   ~    (H+2z)tan6] 


The  boundary  conditions  for  the  additional  terms  alone 
therefore  become: 


x    =  0  at  z  =  +  ? 
rz  2 


a     =  0  at  z  =  ±  | 
z  2 


D0     H 
u  +  u,,  .   ,  .   =  w+w,-,  ,   .  .   =  0  at  r= —  -(■=■+  z)tan( 
flat  plate       flat  plate  z     2 

or 

Dn    h 
u  +  wtane  +  (u+wtan6)^n  ,   .    =0  at  r=—  -  (~   +  z)tan6 

flat  plate         z     2 


For  the  additional  terms  to  satisfy  the  first 
boundary  condition  we  find: 


Trz  =  °  =  (8c  "  8vC  "  6vb  "  144k)  r2j: 


+    (-6b  +  6vb  -  16c  +  8vc  +  144k)  r3 


+  (32g  -  32vg  -  4vf  +  32n)  r 
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Since  the  constants  b,  c,  k,  g,  f,  and  n  are  independent 

3 

of  r,  the  r   and  r  terms  must  each  equal  zero.   Therefore 

144k  =  6b  -  6vb  +  16c  -  8vc 
and 

32n  =  (6b  +  8c)H2  +  (-32g  +  32vg  +  4vf) 


Considering  the  second  boundary  condition  we  require 
that: 

a   =  0  =  (-4b  -  M  C) (±H3)  +  (8q  -  4vq  +  6p  -  6vp) 

Therefore 

6p  (1-V)  =  -8q  +  4vq 
and 

b  =  -  ±c 

D0     H 
At  the  boundary  r  =  -5-  -  (y  +  z)  tan  6 


u  =  ~    [(-16b  +  8vb  -  Mc  +^  Vc)  z4  tan 


+  (6b  -  6vb  +  8c  -  8vc)  z4  tan3( 


+  (8b  -  4vb  +  ^c  -  2y   vc)  Az3 
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+  (-9b  +  9vb  -12c  +  12vc) A2z3tan29  +  (6b  +  8c)H2z2tan( 


+  (|b-|vb+  6c  -  6vc)A2z2tan( 


+    (-32g   +    32vg   -    4f   +    4vf)z2tan9 


+    (-|b+   |vb-    c   +   vc)A3z   +    (-3b   -    4c)H2Az 


+    (16g   -    16vg  +    2f   -    2vf)Az    +    (-2q)    z    tan  9  +    (q)  A  ] 


and 


w  =  ^    [(4b  +   4vb  +  ^-c    +  ^-  vc)    z4 


+  (-12vb  -    16vc)z4tan29    +    (~|b  +   |vb    -    2c   +    2vc)  z4tan49 


+(12vb  +   16vc)Az3tan9    +    (3b   -    3vb  +    4c    -    4vc)Az3tan39 


+  (-3vb-    4vc)A2z2    +    (-fb+   |vb-    3c   +    3vc)  A2z2tan2f 


+(-6b  -  8c)H2z2  +  (4vf  +  32vg)z2  +  (3b  +  4c)H2z2tan2( 


+(16g  -  16vg  +  2f  -  2vf ) z2tan29  +  (-3b  -4c)H2Az  tan 


20 


+    (|.b_  2  Vb'+   c   -   vc)A3z    tan 
4  4 


+    (-16g   +    16vg  -    2f   +    2vf)Az    tan  6  +    (6p   -    12vp  +    8q-8vq)  z 


+  (-Ab+   ^  Vb  _   |c+   |  vc)A4    +    (|b+   c)H2A2 


+     (4g    -    4vc  +   if     -    |  vf)     A2] 


Substituting  for  p  and  b  we  find 


u   =   ij^    [  (-32g   +    32vg  -    4f   +    4vf)z2    tan 


+    (16g   -    16vq  +    2f   -    2vf)Az  (4) 

+     (-2q)     z    tan  0  +    q  A  ] 
and 

1+v  2  ?         2 

w  =  -g-    [  (4vf  +    32vg)  z     +    (16g   -    16vg+2f   -    2vf)z   tan   6 

4v 
+    C-16g   +    16vg   -    2f   +    2vf)Az   tan  9  +   -^  z  (5) 


+    (4g   -    4vg  +   If  -   I  vf)A2] 
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Consequently,  it  is  apparent  that  neither  of  the 
final  boundary  conditions  can  possibly  be  satisfied  for 
the  third  or  fourth  order  terms  using  this  solution  form, 
This  occurs  because  the  first  boundary  condition  causes 
the  fourth  order  constants,  f  and  g,  to  cancel  in  the 
second  boundary  condition.   This  produces  a  relationship 
between  b  and  c  which  prevents  them  from  contributing 
to  the  displacements  u  and  w.   However,  the  addition  of 
higher  order  terms  may  allow  all  boundary  conditions  to 
be  met. 


SOLUTION  WITH  INFINITE  ORDER 
POLYNOMIAL  ASSUMPTION 

In  order  to  investigate  this  possibility  we  consider 
a  stress  function  of  order  n.   Therefore 

,        n  .  ,   n-2  2  ,     n-4  4  . 

<J>   =az  +bz    r  +cz   r  +.... 

Tn    n      n  n 

The  terms  contain  only  even  powers  of  r  since  only  even 

2  2    2 

powers  will    satisfy   equation    (3).      Therefore    V    <j)    and   V    V    <J> 

may   be   expressed   as: 

V2c|>    =    [2-2b     +   n(n-l)a    ]  zn~2    +    [4-4c      +    (n-2)(n-3)b    ]  zn_4r2 
r  n  nJ  n  nJ 


+    [6-6d      +    (n-4)(n-5)c    ] zn-6r4   + 

n  n 


and 


V2V2(J)    =    [2-2-4. 4c   +2-2-2(n-2)  (n-3)b 
T  n  n 


+   n(n-l)  (n-2)  (n-3)a    ]  zn~4   +    [4-4-6-6d  +2- 4- 4  (n-4)  (n-5)cn 


+    (n-2)  (n-3)  (n-4)  (n-5)b    ]zn    6r2    +    [  6  •  6  •  8-  8e  +2  •  6  •  6  (n-6)  (n-7)  d 

n  n  n 


+    (n-4)  (n-5)  (n-6)  (n-7)c   ]  zn    8r4   +    ... 
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2  2 
Consequently,  in  order  to  satisfy  V  V  0  =  0  we  find: 


2(n-2)  (n-3)  .   _  n  (n-1)  (n-2)  (n-3) 
n        4-4       n        2 • 2 • 4 • 4       n 


d     3  (n-2)  (n-3)  (n-4)  (n-5)      2n  (n-1)  (n-2)  (n-3)  (n-4)  (n-5) 
n  4- 4- 6- 6        n  2-2.4.4.6.6  n 


e   =    4  (n-2)  (n-3)  (n-4)  (n-5)  (n-6)  (n-7) 
n  4- 4- 6- 6- 8- 8  n 


3n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)  (n-6)  (n-7) 

2-2-4-4-6-6-8-8  an 


etc. 

o 

Using  Greek  letters  to  distinguish  terras,  V  <j>  =  0  is 

satisfied  if: 

a      _        n(n-l) 


n       2-2     n 


.  n(n-l)  (n-2)  (n-3) 
Tn       2-2- 4- 4        n 


R      =   -n(n-l)  (n-2)  (n-3)  (n-4)  (n-5) 
n  2- 2- 4- 4- 6- 6  n 


etc. 

The  shear  stress  i    at  any  point  can  therefore  be 

rz      2    c 

represented  as: 
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xrz  =    [(-4  +    2v).{(n-2)(h-3)bn.}   +    (-i+v){n(n-1)(^2)(n-3)}an 


+   n(n-l)(n-2)(n-3),  ^]zn-4r   +    [  (  3_v)  {  (n-2)  (n-3)  (n-4)  (n-5)  ^ 


+    to   ^.^  J"'"'1)  (n~2)  ("-3)  (n~4)  (n~5)  \ 
+    (2-v){ _— }an 


n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)  .    n-6    3 
F2M VZ         r 


+    [(-4   +   v)     {  (n"2)  (n"3)  (n"4)  (n"5)  (n"6)  (n"7)  }b 

4' 4*  6  n 


+    (-3   +    V)     ("(n-D  (n-2)  (n-3)  (n-4)  (n-5)  (n-6)  (n-7)  >a 

2  •  2  •  4  •  4  •  6  n 


,    n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)  (n-6)  (n-7)         ,    n-8    5 

2«2-4-4-6  nJZ        r 

etc. 

The   axial   stress    a      can   also   be    represented   at   any 
point   as : 


oz   =    [(2-v)2-2(n-2)b      +    (1-v)  n  (n-1)  (n-2)  a 


-   n(n-l)  (n-2)an]zn    3 


+    [(-3+v)  (n-2)  (n-3)  (n-4)b   +  (-2+y)  n  (n"1  (n~2)  {n~3)  (n"4)    a 

n  2-2  n 


,    n(n-l)  (n-2)  (n-3)  (n-4)  n-5    2 

+  2^2  VZ        r 
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4.  ha    ,o  (n-2)  (n-3)  (n-4)  (n-5)  (n-6) 

/ 

•  fa  ^  n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)  (n-6) 
+  C3-V)  2-2-4-4  an 


n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)  (n-6)      n-7  4 
2-2-4-4  anJZ    r 


etc. 

The  displacements  u  and  w  may  also  be  represented  at 
any  point  as: 


1+v  rr  o  /   o\i_   i  n(n-l)  (n-2)    ,  n-3 

u  =  -=r-  [{-2(n-2)b   +  — i ^ '-   a  }z    r 

t,  n         2        n 


xfn(n-l)  (n-2)  (n-3)  (n-4)     ,  2  (n-2)  (n-3)  (n-4)  . 
+  1        2  ^2^4  an  +  —       4        Dn 


n(n-l)  (n-2)  (n-3)  (n-4)    ,  n-5  3 
2-2-4  ri      r 


,   2n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)  (n-6) 

1  2- 2- 4- 4- 6  n 


3  (n-2)  (n-3)  (n-4)  (n-5)  (n-6) 

4-4-6  n 


+  n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)  (n-6)     ,  n-7  5,      ,,, 

+  ~ — ~ — -j — -. — 7 a  >  z    r  J      (6) 

2« 2-4-4-6  n 


etc. 
and 
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w  =  ±±^[{(l-v)n(n-l)an+(2-v)2-2b  -n(n-l)a    }zn   2 

JZj  .  II  II  II 


+   {  (-2+2V)    n(n-l)(n-2)  (n-3)    a 

z  •  2.  n 


/    -> ,  o    \  /       n\  /       -.M^       ,    n(n-l)  (n-2)  (n-3)  -,    n-4    2 

+    (-3+2v)  (n-2)  (n-3)bn   +   — * ^-^ — ctn}z         r 


,     ,,,    -,»    n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)     n 
+    H3-2V)  2.2-4-4  an 


+    (4-2v)    jn=2)(n-3)(n-4)(n-5)    b 

4.4  n 


n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)  i    n-6    4  ,    . 

2-2-4-4  °VZ        r  u; 

etc. 

These  equations  may  therefore  be  utilized  to  form  a 
series  of  polynomial  solutions  to  equations  (1)  and  (2)  . 
The  ability  of  this  form  to  produce  a  solution  to  the 
desired  equilibrium  and  compatibility  equations  as  well 
as  the  boundary  conditions  may  therefore  be  easily  investi- 
gated for  any  desired  number  of  terms. 

Considering  a  sum  of  general  terms  the  shear  stress 
at  z  =  +-s   becomes: 


Tr-7    =   I    C(-4+2v)  (n-2)  (n-3)b   +  (-l+v)n(n  l]   (n  2)  an 


n(n-l)  (n-2)  (n-3)      f+H%n-4 
2  nJ  K~2} 
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[(3-v)  {n~2)   (n"3)  (n"4)  (n"5)  b 


>  to   „i  n(n-l)  (n-2)  (n-3)  (n-4)  (n-5) 
+  (2-v)  2TYT4 an 


n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)     .  ,+H>n  6r3 
2^"4  V  U2J 


+  r(-4+v)  (n-2)  (n-3)  (n-4)  (n-5)  (n-6)  (n-7)  fa 
L l     '  4 • 4- 6  n 


.       n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)  (n-6)  (n-7)  a 
+  (-3+v)  2-2-4-4-6  an 


-  n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)  (n-6)  (n-7)    ,  ,+H  n  8r5 
2-2-4-4-6  nJ  K~2! 


etc. 

In  order  to  satisfy  the  first  boundary  condition,  t    =0 
and  therefore  each  order  of  r  must  be  independent  and  must 
equal  zero.   It  can  also  be  seen  that  even  values  of  n  do 
not  change  sign  for  z  =  ±-~   but  odd  values  of  n  do  change 
sign.   As  a  result,  for  each  order  of  r  the  even  and  odd 
terms  must  be  independent  and  must  each  equal  zero. 
Expressing  each  sum  in  terms  of  its  lowest  order  member 
we  find: 


For  r 
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a4   "11=4,6,8. 


2(-4+2v)  (n-2)  (n-3)        ,H.n    4    , 
(1-v)    4-3-2-1'  K2'  n 


+  E 


11=4,6,8, 


n(n-l)  (n-2)  (n-3)        ,H.n    4 
(1-v)    4-3-2-1  V  n 


(8) 


+  £ 


n=6,8,10 


-n(n-l)  (n-2)  (n-3)     ,H   n"6 

4-3-2-1  V  n 


_y  2(-4+2v)  (n-2)  (n-3)       H ,n~5    . 

a5    ~  11=5,7,9...        (1-v)     5-4-3-2  V  n 


+  Z 


n(n-l)  (n-2)  (n-3)     (H}n   5 

2  n 


n=5,7,9...     (1-v)       5-4-3-2 


(9) 


+ 


n=7,9,ll. . 


-n(n-l)  (n-2)  (n-3)     ,H  n~7 

5-4-3-2  v2;  n 


For    r" 


1      =  E 
6         n=6,8,10. 


-2-2  (3-v)  (n-2)  (n-3)  (n-4)  (n-5)     ,H,n-6 

(2-v)     6-5-4-3-2-1  (2}  n 


+  1 


n=6,8,10. 


n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)     (H.n    6  (       . 

(2-v)     6-5-4-3-2-1  2  n 


+  E 


n=8, 10,12. . 


-n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)     ,H    n    8 

6-5-4-3-2-1  (2'  n 
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l7        n=7,9,ll.. 


•2-2  (3-v)  (n-2)  (n-3)  (n-4)  (n-5)     ,H.n    7 

(2-v)     7-6-5-4-3-2  K2}  n 


j  n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)     ,H 

n=7,9,ll...  (2-v)     7-6-5-4-3-2  ' 21 


n-7 


an         (11) 


+  - 


n=9,ll,13. . 


-n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)     ,H.n   9 

7- 6* 5- 4- 3-2  K2'  n 


For   r" 


n-i 


=  v  2-2(-4+V)  (n-2)  (n-3)  (n-4)  (n-5)  (n-6)  (n-7)    H> 

8         n=8,10,12...  (3-v)     8 • 7- 6  -  5- 4 • 3- 2 • 1  K2}  n 


n-i 


+  S  n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)  (n-6)  (n-7)  .H 

n=8, 10,12..  (3-v)     8-7-6. 5-4-3- 2-1  V  n    l       ; 


z  -n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)  (n-6)  (n-7)    H  n    10 

n=10,12,14..  8-7-6-5-4-3-2-1  K2>  n 


a9    n=9,ll,13, 


2-2(-4+v)(n-2)  (n-3)  (n-4)  (n-5)  (n-6)  (n-7)   ,H   n~9 

(3-v)     9- 8- 7- 6-5- 4- 3- 2  ^2"'  n 


n=9,ll,13. 


n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)  (n-6)  (n-7)     ,H    n    9  ,       . 

(3-v)     9-8-7.6-5.4- 3-2  (2~'  n    K±J> 


n=ll,13,15 


-n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)  (n-6)  (n-7)  ,H   n_11 

9-8-7-6-5-4-3-2  V  an 


etc, 
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When  higher  order'  equivalents  of  the  a  terms  are 
substituted  into  the  lower  order  equivalents  we  find  that 
each  a   can  be  expressed  as  a  sum  of  equal  and  higher 
order  terms  of  b  and  a.   In  the  event  that  an  a  is  zero, 
the  same  order  b  can  be  expressed  as  a  sum  of  an  equal 
order  term  of  a  and  higher  order  terms  of  b  and  a. 
Similarly,  if  an  a  and  b  were  both  zero,  the  same  order 
term  of  a  can  be  expressed  as  a  sum  of  higher  order  terms 
of  b  and  a. 

Disregarding  the  flat  plate  solution  and  attempting 
a  general  solution,  the  boundary  conditions  for  a   become: 

a   =  0  at  z  =  Tt- 
z  2 

a   =  -P  at  z  =  -77 
z  2 

2    4 
Consequently,  the  sums  of  a      terms  containing  r  ,  r   etc. 

must  all  equal  zero  at  z  =  i-r  whereas  the  sums  of  terms 

which  are  functions  only  of  z  must  be  non-zero.   The  odd 

H  H 

order  terms  do  not  change  sign  for  z  =  -~    or  z  =  — =-,  however, 

the  even  order  terms  do.   The  odd  and  even  sums  containing 

2    4  H 

r  ,  r   etc.  therefore  each  equal  zero  at  z  =  t-j.      The  sum 

of  odd  terms  which  is  only  a  function  of  z  must  therefore 

p 
be  equal  to  -■=■  while  the  sum  of  even  terms  must  be  equal 

P         HP  H 

to  -j  at  z  =  —  and  — *-  at  z  =  -~-.   This  results  in  the 

following  equations : 
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r  [(2-v)2-2(n-2)b;  +  (l-v)n  (n-1)  (n-2)  a 

1   C  1  n 

n=3, 5, 7 


h  n~3    P 
-  n(n-l)  (n-2)  an]-  (|)     =  ~ 


y  [(2-v)2-2(n-2)b  +  (l-v)n (n-1)  (n-2) an 

n=4, 6,8.. 


-  n(n-l)  (n-2)  aj  (§)     =  | 


2 

For  r 


j         [  (-3+V)  (n-2)  (n-3)  (n-4)bn  +  .(-2+v)n  (n- 1)^2^-3)^-4)^ 
n=5, 7,9. . 


,  n(n-l)  (n-2)  (n-3)  (n-4)     ,H.n  5    n 
2-2  nJ  l2J 


E     [  (-3+v)  (n-2)  (n-3)  (n-4)  b+  (-2+v)  "  <£="   (n-2)  (n-3)  (n-4) 
n=6,8,10... 


+  n(n-l)  (n-2)  (n-3)  (n-4)    .  ,H.n  5 
2-2  nJ  v2; 


etc. 

If  the  values  of  "a"  satisfying  the  shear  stress  boundary 

condition  are  substituted  we  find  for  a.,: 
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-1  (-2+v) 2  .      1 

a3  "  2(l-v)  3-2-1  P  +  (1-v) 3    3  +  1-V   a3 


2»3  [-2(-2+v)]   ,H  2 
(1-v)  3-2-1     V    5 


y  ,2 (n-2)  [-2(n-4)  (-2+v)]  ,H  2 

n=7,9,ll. .{  (1-v)  3-2-1        V 


-2(3-v)  (n-2)  (n-3)  (n-4)  (n-5)  }     ,H  n~5 

(2-V)  3-2-1-2-3  v2;      n 


,  i  n(n-l)  (n-2)  [2-(n-3)l  ,H> 

n=7,9,ll..1   (1-v)  3-2-1-2      y2l 


,  n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)  ,  ,H  "  5 

(2-v)  3-2-1-2. 3-2      iK2>  n 


+  E  -n(n-l)  (n-2)  (n-3)  (n-4)  (n-5)  ,H  n~7        ,   . 

n=9,ll,13..        3-2-1-2-3-2  v2;      n     v   ; 


etc. 

By  substituting  the  appropriate  orders  of  a,  a_.  may 
be  expressed  in  terms  of  P,  b_.,  a_,  b  ,  and  higher  odd 
orders  of  b  and  a.   The  a   term  does  not  appear  directly 
in  this  expression  although  all  orders  of  a  do  appear. 

Similarly,  we  find  an  expression  for  fourth  and  higher 
orders  of  the  r-independent  sum: 
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2 
(1-v)  6-5-4  [(|) '  -1]  (|)  a6 

n-3 

+  y        2(n-2)  [2(2-v)  +  (-4+2v).(n-3)  ]  (§)     b 
n=6,8,10... 

n-3 
+E  (-n)  (n-1)  (n-2)  [(n-3)-l]  (§)     a 

n=6,8,10...  ^ 

+  Z        (l-v)n(n-l)  (n-2)  t(|)   -1]  (f)     a   =|    (15) 
n=8,10,12...  z         z      n 


In  this  boundary  condition  the  shear  stress  boundary 
condition  causes  all  fourth  order  terms  to  cancel  resulting 

in  an  expression  for  sixth  or  higher  order  terms. 

2 

The  odd  sum  for  terms  which  are  functions  of  r 

becomes : 


h  =1  -n(n-l)  (n-2)  (n-3)  r      (1-v)  (n-4)  +  (-2+v)     .  H  n  5 

5n=5,7,9...   2-2-3      L (-3+v) 2 (1-v)  +  (-2+v)  (-4  +  2v) J  { 2'  n 


-(n-2)  (n-3)  ,  (-3+v)  2  (1-v)  (n-4)  +  (-2+v)  (-4  +  2v)  .  ,H  n  5 
+n=7  9  11  3'2       (-3+v)2(l-v)  +  (-2+v)  (-4+2v)    M2>    Dn 

(2-v)n(n-l)  (n-2)  (n-3)  (1-v)      (r_4)  g,2^  H,n"7a 


i     2-2- 3[  (-3+v)  2(l-v)  +  (-2  +  v)  (-4  +  2v)  ]  L  v    '  K  2'       iJV2'     n 
n=7,9,ll 

(16) 
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It   is  apparent  by  considering  the  higher  order  terms  that 

none  of  these  terms  become  zero.   The  even  sum  for  terms 

2 
which  are  functions  of  r   becomes: 


n-8 

a 


n 


b  =  e  n(n-l)  (n-2)  (n-3)  (n-4) 

8    n=8,10,12...     6-5-4-3-2-2 


r3(3-v)  (n-6)+4(-2+v)  (n-5)  (n-6)  (n-7)  ,  ,H. 
1    (-2+v)  (-2)  (-4+v)+2(3-v)  (3-v)    J  K2] 


+E  -(n-2)  (n-3)  (n-4)  (n-6) 

n=10,12,14. . .     6-5-4-3.2-1 


r(3-v)23-2>l+(-2+v)  (-2)  (-4+v)  (n-5)  (n-7)  .  ,H.n~8 

C-2+v)  (-2)  (-4+v)+2(3-v)  (3-v)      J  K2]  n 


-(3-y)  (-2  +  y)n(n-l)  (n-2)  (n-3)  (n-4)   n  ,  _^wH.n"8 
E    6.  5- 4.  2- 2  [(-2  +  v)  (-2)  (-4+v)  +2  (  3-v)  2]  U~ln  b;  J  [2}         an 
n=10,12,14. . . 


(17) 


Once  again  the  lowest  order  terms  cancel.   The  same  type 


behavior  is  found  for  terms  which  are  functions  of  r4,  r6 , 

8 
r  ,  etc. 

Several  important  characteristics  appear  in  these 
equations.   The  odd  order  terms  of  fifth  order  or  higher 
result  in  an  expression  of  b  in  terms  of  the  same  and 
higher  orders  of  a  only.   Since  the  third  order  terms  did 
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not  appear  in  the  shear  stress  boundary  condition,  a~  is 
expressed  in  terms  of  P,  b.,,  and  a  terms  of  third  order  and 
higher. 

The  even  terms  behave  similarly  to  each  other  for 
each  order  of  .?.   In  each  case,  the  lowest  order  terms 
cancel  each  other  and  the  next  higher  order  of  b  can  then 
be  expressed  in  terms  of  the  same  and  higher  orders  of  a. 

These  derived  expressions  for  the  b  terms  can  be 
substituted  into  the  expressions  from  the  first  boundary 
condition.   The  final  result  for  third  order  will  be  an 
expression  for  a_  in  terms  of  b_,  P,  and  third  and  higher 
orders  of  a.   For  fifth  order  and  higher,  the  result  will 

be  expressions  for  a   and  b   in  terms  of  a   of  order  n  and 

c  n      n 

higher.   The  final  result  for  fourth  order  terms  is  an 

expression  for  a.  in  terms  of  b  and  a  of  fourth  order 

and  higher  since  the  fourth  order  terms  cancelled  out  of 

the  final  boundary  condition.   For  sixth  order  terms,  the 

final  result  is  expressions  for  a^  and  b^  in  terms  of  P 

6       6 

and  a  of  order  six  and  higher.   For  eighth  and  higher  order 

even  terms,  the  final  result  is  expressions  for  a   and  b 

r  n       n 

in  terms  of  a  of  order  n  and  higher.   Consequently,  the 

unevaluated  constants  available  to  satisfy  the  final 

boundary  condition  are: 

b_   a0   b.   a .   ar   a,   a-   a0  . . . . 
3,        3 ,   4,   4,   5,   6,   7,   8 
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The  above-mentioned  behavior  of  the  equations 
satisfying  the  first  boundary  conditions  combined  with  the 
previous  attempted  solution  using  third,  fourth,  and 
sixth  order  terms  allows  determination  of  the  crucial 
factors  which  will  or  will  not  permit  satisfaction  of  the 
final  boundary  condition.   In  the  previous  attempted 

solution  equivalents  of  b..,  a^,  b . ,  and  a.  appeared  in 

2 

the  z  ,  z,  and  constant  terms  of  u  and  w.   However,  each 

term  contained  the  same  relationship  between  b   and  a_  and 
between  b.  and  a  .   As  a  result,  only  2  of  these  constants 
are  useful  in  satisfying  the  final  boundary  condition.   The 
cancellation  of  the  sixth  order  terms  in  u  and  w  combines 
with  this  to  produce  a  total  of  five  equations  requiring 
satisfaction  but  only  two  constants  with  which  to  satisfy 
them.   As  a  result,  three  additional  constants  are  needed 
in  order  to  satisfy  the  final  boundary  condition. 

From  equations  (6),  (7)  ,  (9)  ,  and  (16)  it  can  be  seen 
that  the  addition  of  fifth  order  terms  will  introduce 
another  constant  a   to  assist  in  satisfying  the  final 
boundary  condition.   In  order  to  satisfy  the  final 
boundary  condition,  therefore,  two  higher  order  constants 
are  needed. 

Since  the  sixth  order  terms  appear  in  equations 

(6)  and  (7)  as  a.,  and  P  and  since  the  a,  coefficients 

d  6 

cancel,  higher  order  terms  are  needed  to  cancel  the  P 
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o       4 
multiples  of  z   and  z  .   The  next  available  terms  are 

seventh  and  eighth  order  but  equations  (6)  and  (7)  show 

5 
that  their  introduction  will  cause  u  and  w  to  have  z 

and  z   components.   The  z   term  will  contain  only  aR 

5 
and  the  z   term  wxll  contain  only  a_  and  afi. 

If  a_  and  a0  are  utilized  to  eliminate  the  P 

3       4        '5       6 
multiples  of  z   and  z  ,  the  z~  and  z   terms  will  be 

unsatisfied.   The  introduction  of  higher  order  terms  will 

not  solve  this  problem  since  equations  (8)  through  (17) 

show  that  all  higher  order  terms  behave  as  the  seventh 

and  eighth  order  terms  and  the  result  would  just  be  the 

raising  of  the  order  of  the  unsatisfied  terms  in  the  final 

boundary  condition.   The  only  possible  way  to  completely 

satisfy  the  final  boundary  condition  would  require  that 

the  coefficients  of  two  of  the  higher  order  terms  cancel 

as  the  coefficients  of  ar    cancelled. 

b 

Investigating  this  possibility  for  a     we  find  for 

o 

the  z   term  in  u: 


u=  1+V  t_  8_^_6  3(2_v)      z6  tan( 


By  examination  of  equations  (8)  through  (17)  it  can  be 

seen  that  none  of  the  coefficients  of  a   will  cancel  in  u 

n 

or  w  for  n  ^  7.   As  a  result  it  is  impossible  to  satisfy 
all  the  boundary  conditions  for  this  problem  with  even  an 
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infinite  series  of  polynomials. 

Since  the  possible  values  of  z  are  all  fairly  small 
compared  to  the  radius,  the  components  of  u  and  w  at  the 
boundary  from  the  two  unsatisfied  highest  orders  should 
approach  zero  as  the  order  increases.   Consequently  it 
would  seem  that  by  taking  a  sufficiently  high  order 
equation  and  satisfying  the  boundary  conditions  for  all 
but  a  negligible  value  of  u  and  w,  an  accurate  approximation 
can  be  obtained. 

This  approach  produces  a  sufficiently  small  error  in 
the  final  boundary  condition  but  when  the  resulting  con- 
stants are  used  to  calculate  the  deflection  of  the  center 
of  the  viewport,  it  does  not  compare  well  to  the  experi- 
mental results  obtained  by  Stachiw.   This  occurs  because 
although  only  the  high  order  boundary  condition  is 
unsatisfied,  the  high  order  term  affects  the  evaluation 
of  all  lower  order  constants.   This  is  apparent  from 
equation  (14).   Consequently,  this  method  cannot  even 
provide  a  reasonable  approximation. 


DISCUSSION  OF  RESULTS 

A  solution  form  for  a  viewport  in  a  rigid  mount  is 
proposed  which  appears  to  have  an  ample  number  of  terms 
to  be  able  to  satisfy  all  boundary  conditions.   The 
compatibility  equations  are  such  that  the  satisfaction  of 
the  first  two  boundary  conditions  couples  the  available 
terms  so  that  the  final  boundary  condition  cannot  be 
satisfied.   Although  for  high  orders  the  error  is  small, 
the  terms  are  coupled  in  such  a  way  that  all  terms  are 
affected. 

The  inability  of  this  approach  to  provide  satisfactory 
results  is  due  to  the  compatibility  conditions  and  by 
the  assumption  of  plane  stress.   The  nature  of  the 
equations  with  plane  stress  is  such  that  it  cannot 
accommodate  a  slanting  surface  for  the  final  boundary 
condition. 

Consequently,  even  for  thin  viewports  plane  stress 
cannot  be  assumed  in  forming  an  analytic  solution.   The 
general  three-dimensional  compatibility  equations  must  be 
utilized  and  the  simplification  of  stress  functions  cannot 
be  used. 

The  equilibrium  and  compatibility  equations  can  best 
be  solved  by  expressing  them  in  terms  of  u  and  w.   The 
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equilibrium  equations  become: 


and 


(G+A)  si  (I?  +  ?)  +  G(T"2  +  F  ^}  +  t^v  G  "X)  TT  =  ° 

or  ■  &z 


The  compatability  equations  become 


A  (1-v)  {&3u   +  l^u+   6~3u 
v     <sr3    r  6r2    6r6z2 


+  a(1^__L1h  +  JL+   63w   +  1  J^w.   1  S2u    63w 
r  6r2    r2  6r    r3    6r26z  '  r  6r6z    r  ^2    ^3' 


_  4G   6u  _  u     X  ,1  62u  _  2  6u    2u  ,  6 3u  .   6 3w     . 
r2  ^    ^  +  V  (r  ^2    ?  ^   ?   ^3  +  ^  =  ° 


1-v  ,1  ^fu    _1_  6u  ,  _u_    1  62u. 
V    r  6r2  "  r2  6r    r3    r  ^2> 


SrJ  Sr  6r6z         6r25z        r    ^6z         6z3; 


+  i§    (|H  _  H)    +  i    (1  £!u  +  J_  ^u  _  _u_       l    52w   )    _    . 
2     v6r        r'  v     vr    r    2    +      2    A-r  5    +   7  TTT^)    "    ° 


r 


^      r    6r2         r2    fir        r3         r    6r6: 
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1-v  ,  *3w  1  52w    s\ 

v   ^6r26z    r  6r6z    6z3J 


r  6r2    r2  6r  '  r3    sr3    r  Mz    M  2* 


V   r  6z2    6r6z2    6Z3- 


G  T2  -!!«-  +  ¥*   +  1  62u     2  Su 
6r26z    6z3  +  rfi^-l^ 


+  6\  +  2  _5\     i  fiw  =  Q 
6r      5r6z2   r2  6r 


Once  a  solution  form  is  found  which  satisfies  these 
equations,  the  boundary  conditions  can  be  utilized  to 
evaluate  the  constants  in  the  equation. 


RE  COMMEND AT I  ON  S 

The  solution  of  the  aforementioned  equilibrium  and 
compatibility  equations  will  not  be  simple.   A  pre- 
liminary investigation  reveals  that  no  simple  assumptions 
provide  solutions.   If  no  analytic  solution  is  obtained, 
designs  can  still  be  safely  made  using  numerical 
techniques  but  an  analytic  solution  would  prove  much 
more  useful  to  others.   An  analytic  solution  would  also 
provide  a  great  deal  of  insight  into  the  stresses  in  the 
reinforced  area  of  the  hull  which  serves  as  a  mounting. 
For  these  reasons,  further  work  should  be  done  to  try  to 
achieve  a  solution  to  this  problem. 
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